Equations of motion and conservation laws in a theory of stably stratified turbulence 
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The letter considers non-isothermal fluid flows and revises simplifications of basic hydrodynamic 
equations for such flows arriving eventually to a generalization of the Oberbeck-Boussinesq approx- 
imation valid for arbitrary equation of state including both non-ideal gases as well as liquids. The 
proposed approach is based on a suggested general definition of potential temperature. Special at- 
tention is put on the energy conservation principle, and it is shown that the proposed approximation 
exactly preserves the total mechanical energy by approximate equations of motion. The principal 
importance for any turbulent boundary layer model to respect the conservation laws is emphasized 
explicitly. 



Introduction 

In this letter, which has partially a pedagogical character, 
we reconsider simplifications of the basic hydrodynamic 
equations of motion for non-isothermal fluid flows, which 
can serve as a basis for a consistent theory of atmospheric 
turbulent boundary layers with stable temperature strat- 
ification. The reason for doing so stems from the fact 
that textbooks derivations of the celebrated Oberbeck- 
Boussinesq approximation are too specialized either to 
fluids which density variation is very small or to ideal 
gases and cannot be applied, for example, to humid air. 
Our generalization of the Oberbeck-Boussinesq approx- 
imation is based on a suggested below generalized def- 
inition of the potential temperature, which is valid for 
arbitrary equation of state, while its standard definition, 
widely used in meteorology, is applicable only for ideal 
gases. 

We discuss also a form of mechanical energy (consist- 
ing of the kinetic and potential energies) that is exactly 
preserved by the resulting approximate equations of 
motion. This integral of motion plays a crucial role 
in modeling stably stratified atmospheric turbulent 
boundary layers (see our paper in the same issue), 
allowing to control further approximations with the goal 
of a reasonable, simple, but still adequate description 
of turbulent fluxes that determine basic physics of the 
turbulent boundary layers. 



1. Generalization of the Oberbeck-Boussinesq 
approximation for non-ideal gases and liquids 



1A. Basic hydrodynamic equations. The system of 
hydrodynamic equations describing a fluid in which the 
temperature is not uniform consists of the Navier-Stokes 
equation for the fluid velocity, U(r,t), a continuity equa- 
tion for the space and time dependent (total) density of 
the fluid, p(r, t), and of the heat balance equation for the 



(total) entropy per unit mass, S(r,t), [l|: 



dt 



p(—+U-v)u = -Vp + g p+V ■ pVU 

V-( P U) = 0, 
U-v)s = V • kVS . 



dp 
at 

fB 



(1.1a) 
(1.1b) 
(1.1c) 



Here p is the pressure, g = —zg is the vertical accelera- 
tion due to gravity, p and k are the (molecular) dynam- 
ical viscosity and heat conductivity. 

These equations are considered with boundary condi- 
tions that maintain the solution far from the equilibrium 
state, at which U = S = 0. These boundary conditions 
are U = at zero elevation, U — const at a high eleva- 
tion of a few kilometers. This reflects the existence of a 
wind at high elevation, but we do not attempt to model 
the physical origin of this wind in any detail. The only 
important condition with regards to this wind is that it 
maintains a momentum flux towards the ground that is 
prescribed as a function of the elevation. Similarly, we 
assume that a stable temperature stratification is main- 
tained such that the heat flux towards the ground is pre- 
scribed as well. In the entropy balance Eq. lc|) we have 
already neglected the viscous entropy production term, 
oc p\VU\ 2 , assuming that the temperature gradients are 
large enough such that the thermal entropy production 
term on the RHS of Eq. p.lc[) dominates. Actually, this 
assumption is very realistic in meteorological and oceano- 
graphic applications. For simplicity of the presentation 
we restrict ourselves by relatively small elevations and 
disregard in Eq. (|l.la[) the Coriolis force (for more de- 
tails, see [2j). 

On the other hand we assume that the temperature 
and density gradients in the entire turbulent boundary 
layer are sufficiently small to allow employment of local 
thermodynamic equilibrium. In other words, we assume 
the validity of the equation of state, and that the entropy 
S is a state function of the local values of the density and 
pressure: 



P = P{ T ,P) , S = S(p,p) . 



(1.2) 
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In the same manner we will neglect the temperature de- 
pendence of the dissipation parameters v and k. 
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Pressure fluctuations caused by turbulent velocity 
fluctuations u propagate in a compressible medium 
with the sound velocity c s , causing time dependent 
density fluctuations of the order of (u/c s ) 2 po, where 
Pa is the mean density. Assuming that the square of 
the turbulent Mach number M% = (u/c s ) 2 is small 
compared to unity, we can neglect in Eq. (|l.lb[) the 
partial time derivative, hence V • {pU) — (see e.g. 
Pi). Even tropical hurricanes of category five have the 
mean wind velocity U is below 300 Km/h. Usually, the 
turbulent velocity fluctuations u are less then U/10, 
i.e. even in these extreme conditions u < 30 Km/h and 
Ml < 10- 3 (with c s ~ 1200 Km/h). Therefore the 
incompressibility approximation V • {pU) = is well 
justified in atmospheric physics. In the ocean where the 
sound velocity is even larger and water velocities even 
smaller, this approximation is quite excellent. 

IB. Isentropic basic reference state. In quite air, 
without turbulence, the pressure and the density depend 
on the elevation z simply due to gravity. For example, 
in full thermodynamic equilibrium the temperature is 
uniform, z-independent, and the density decreases expo- 
nentially with the elevation. However, this equilibrium 
model of the atmosphere is not realistic, and cannot be 
used as a reference state about which the actual dynam- 
ics is considered. A much better reference suggestion is 
a state in which the entropy is space homogeneous. In 
this model the thermal conductivity (leading to the tem- 
perature homogeneity) is neglected with respect to heat 
transfer due to the vertical adiabatic mixing of air, lead- 
ing to a z-independent entropy. We refer to the isentropic 
model as a "basic reference state" and denote this state 
of the system with a subscript " b" : 

5 b = S(p h ,p h ) = const, p h = p(T b ,p b ) ■ (1.3) 

The first of Eq. (|1.3[) relates the gradients of the pressure 
and density in this state: 



= VS b 



dpb 



dpb 



(1.4a) 



Another relation between pb and pb follows from the con- 
dition of hydrostatic equilibrium: 



(1.4b) 



Equations (|1.4[) together with the first of Eqs. (|1.3p 
determine the density, pressure and temperature profiles 
in the isentropic basic reference state. 

1C. Hydro dynamic equations in generalized 
Oberbeck Boussinesq approximation 

a. Equations of motion. Denote the deviations of 
the total density, pressure, temperature and entropy from 
the basic reference state as follows: 



P = p-Pb, p = p-pb, 
T = T — Tb , S = S- S h ■ 



(1.5) 



Following Oberbeck 0] and Boussinesq Q, assume that 
these deviations are small: p <S Pb , P Pb- Then one 
simplifies the full system of hydrodynamic Eqs. (jl.ip and 
rewrites them in terms of the fluid velocity 1A and the 
small deviations p and S (instead of p) . The first step is 
very simple: because of Eq. (jl.4b|) 



- gp = Vp - g p . 



(1.6a) 



Next we should relate the deviations p, p and S. In the 
linear approximation Eq. (jl.2[) yields: 

(1.6b) 



S = (dSb/dpb) p p+ (dS b /dp h ) p p 
With the help of Eqs. (|1.4p this gives: 

g p = p-p — T h s . 

Pb c p 



(1.6c) 



Here (3 = g (3 is the buoyancy parameter, f3 is the thermal 
expansion coefficient and c p is the isobaric specific heat 
in the basic reference state (for thermodynamic relations 
see, for example, |5j): 



1 

Pb 



&= iSr 



dpb 
dT h 



PbTb 



dpb 
dSb 



7^ • (l-6d) 



Now Eq. (|1.6ap in the linear approximation yields: 



Vp - gp = Pb 



c. 



, Pb / 

Then Eqs. p. lap can be approximated as 



(1.6e) 



(3— S+— V-pVU . (1.7) 

Co Pb 



K Ph / >-p 

b. Generalized potential temperature. To proceed, 
we generalize the notion of potential temperature O (see, 
e -g- @]) which is traditionally defined as the temperature 
that a volume of dry air at a pressure p{z) and tempera- 
ture T{z) would attain when adiabatically compressed to 
the pressure p„ that exists at zero elevation z — 0. This 
potential temperature can be explicitly computed for an 
ideal gas with the result 



e(z) = T„(p„/ P (*)) 



(7-l)/7 



(1.8) 



where 7 is the ratio of isobaric to isochoric specific heats, 
7 = Cp/c v , and T* is the temperature at zero elevation. 

We want to generalize the notion of the potential tem- 
perature for an arbitrary stratified fluid requiring that in 
the isenotropic basic reference state it would be constant 
O* = T*. A second requirement is that the definition 
will agree with Eq. (|1.8p for an ideal gas. Accordingly 
we define 

e(z) = T*e X p[(S(z)-S h )/c p ] . (1.9) 

For more details see also Q- Indeed, if we employ the 
equation of state and the equation for the entropy of an 
ideal gas, i.e. 

p = pT, S = In (p Cv /p Cp ) + const , (1.10) 

one can easily check that Eq. (|1.8p is recaptured. 
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c. Resulting equations. For small deviations of 
from the basic reference state value T», i.e. up to the 
linear order, Eq. (| 1 .9(1 gives: 

6 d = 6-T, =T»S/c p . (1.11) 

Now we can present Eqs. (jl.lal) . (|l.lb|) (with dp t /dt = 0, 
as explained) and <|1. lc|) as follows: 

t^+U-V^U = -v(^j-/3Q d + vAU, (1.12a) 
V-(a,W)=0, (1.12b) 
+ W v)e d = X A9 d . (1.12c) 

The dissipative terms are important only in the narrow 
region of the viscous sublayer, where we can safely ne- 
glect the z-dependence of ph, /i and k, and consider the 
dynamical viscosity v = p/pb and dynamical thermal 
conductivity x — K l Pb as some z-independent constants. 

Note that for turbulence in liquids (water, etc.) one 
can simplify these equations further. There one can ne- 
glect the effect of adiabatic cooling (together with the 
compressibility), and simply use another reference state 
with constant temperature and density: 

T = T„, pb = P*, Ph = p„ + gp^z . (1-13) 

For this reference state, the standard reasoning (see, e.g. 
[l|) yields the same equations as Eqs. (|1.12[) in which 
again /3, is given by Eq. (|1.6d[) and it is a parameter 
characterizing a particular fluid. In this case ph = p*, 
independent of z, and the potential temperature = T, 
such that d = Td is a deviation of the total temperature 
T from its ground (bottom, or whatever) level T*. 

In the suggested Eq. (|1.12|) the situation is more 
general since we do not assume that the reference state 
has the simple form (| 1 . 1 3[) (with pb = const). Impor- 
tantly, on the RHS of Eq. (|1 . 12a|) the density Pb(z) is 
operated on by the gradient, and the buoyancy term 
— /30 d involves d ^ Td, the deviation of the potential 
temperature defined by Eq. (jl.lip . This definition 
for liquids has nothing in common with the standard 
meteorological definition (|1.8[) . Notice also that for an 
ideal gas (3 = 1/T, and Eq. (|1.12a[) coincides with that 
suggested in the book Q. 

2. Conservation of total mechanical energy 

2A. Dynamical integral of motion. It is important to 
realize that the approximate Eqs. (|1.12[) exactly conserve 
an approximate expression for the total mechanic energy 
of the system in the dissipation- less limit. This total 
energy is the sum of the kinetic, £ K , and the potential 
energy £ P (calculated in the basic reference state): 

£ K = fdrpj^j-, £ P = I dr Pb /3rQ d . (1.14) 



One can check by direct substitution that this sum of en- 
ergies is conserved by equations of motion (Eqs. (|1.12[) ) 
when v — x — 0, i.e the sum £ K + £ P is a dynamical 
integral of motion. 

2B. Statistical integral of motion. We show now 
that the potential energy of a stratified turbulent flow 
£ P , Eq. (|1.14|) . can be presented in a turbulent regime 
as a sum of (time-independent) potential energy of the 
basic reference state, £ P and a "turbulent" potential en- 
ergy, associated with temperature fluctuations, £ P , as it 
was first mentioned by L. F. Richardson Q. Actually, it 
is more instructive to discuss this issue in a more general 
case, when the stratification is caused by some "internal" 
parameter of the fluid, £, not necessarily the potential 
temperature. It can be salinity of water in a sea, humid- 
ity of air, a concentration of particles co-moving with the 
fluid as Lagrangian tracers, etc. 

In the general case then the equation for the potential 
energy of a stratified fluid has the form: 

£ P = g J p(r)z dx dy dz , (1.15a) 

In the basic reference state, the potential energy reaches 
its minimum value referred to as the basic potential en- 
ergy, £ P : 

£ P = g J p b (z)z dx dy dz . (1.15b) 

Clearly, in the equilibrium the fluid density p b (z) de- 
creases with the elevation: dp h /dz < 0. In a turbu- 
lent state the density deviates from its reference value: 
p(r,t) = p b (z) + p(r,t) and the mean potential energy 
(£ P ) exceeds £ P : 

(£p)=£p+£ P . (1.15c) 

We compute the turbulent potential energy, £ P , in the 
case when the internal parameter £ (temperature, etc.) is 
co-moving with the fluid element as a Lagrangian marker. 
In the Lagrangian approach we can consider p as the La- 
grangian marker and introduce a variable z(p,t), which 
is understood as an elevation of the fluid element with 
the density p. Noticing that zdz — ^dz 2 , and integrating 
Eq. (|1.15a[) by parts with respect of z 2 , we can present 
(£ P ) in the Lagrangian approach as: 

<£p> = -§ / ([z(p,t)] 2 )dxdydp . (1.16) 

As a result of turbulent motion, the elevation z(p,t) at 
given p fluctuates and can be decomposed into the mean 
and fluctuating parts: 

z( P> t)=z b (p)+z(p,t), (z(p,t)) = 0. (1.17) 

The substitution of z 2 (p, t) = z 2 + 2z b z + z 2 in Eq. (fTTTTJT) 
leads to three contributions to the potential energy. The 
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first one (originating from z 2 ) describes the basic poten- 
tial energy, Eq. (|1.15b[) . The second contribution, which 
is linear in z, disappears because (z) = 0. The last one 
describes the turbulent potential energy: 



| / ([z(p,t)] 2 ) dxdydp 



(1.18) 



Relating the density fluctuations, p, around the basic ref- 
erence state density profile p b (z) with z 



6p ^^lsz. 
dz 



(1.19) 



and returning back to the Eulerian description in 
Eq. (|1.18|) . one has: 



■f / ^ 



dp b (z) 



dz 



dx dy dz 



(1.20) 



Here we used the transformation formula, similar to 
Eq. (fTTT9| : dp = [dp h (z)/dz] dz. Equation (fj~20f allows 
one to introduce a local density of turbulent potential 
energy per unit mass, 



dp h (z) 



dz 



P 



(1.21a) 



(1.21b) 



such that 

£ p = J p b E P dx dy dz . 
For the particular case of temperature stratification 

= gp, (1.22) 



dp h (z) _ ~dQ 
dz ^ b dz 



P = Ph 



where O — (Od) is the mean potential temperature, 
9 = 9d — O is the fluctuation of potential temperature, 
and the turbulent potential energy per unit mass for the 
stratified turbulent boundary layers is: 



E P — (3Eg/Se , 
where Eg = {9 2 )/2 and S e = dB/dz. 



(1.23) 



Clearly, a consistent statistical description of turbulent 
flows must conserve (in the dissipationless limit) the to- 
tal mechanical energy, which for stratified flows consists 
of kinetic energy of the mean flow, kinetic energy of tur- 
bulent velocity fluctuations, u, and turbulent potential 
energy, E P . To respect this conservation law one has to 
discuss explicitly the balance equation for the turbulent 
kinetic energy that includes components of the Reynolds 
stress tensor (uiUj), and the balance equation for the 
potential energy, which is proportional to (9 2 ). Unavoid- 
ably, the turbulent heat flux vector (u6) is involved into 
the game and requires a separate balance equation for 
itself. Therefore, a consistent statistical description of 
the stratified turbulent flows demands an explicit con- 
sideration of the whole set of the second-order one-point, 
simultaneous (cross)-correlation functions: (iiiUj), (u9) 
and (9 2 ). Our version of such a consistent treatment is 
presented elsewhere in this issue. 

Summary. We presented here a generalization of 
the Oberbek-Boussinesq approximation for temperature 
stratified flows. We showed via a detailed discussion that 
this approach is acceptable for fluids where the local ther- 
modynamic equilibrium is established i.e. the equation 
of state is valid. Thus the applicability of the approxima- 
tion is much wider than for just ideal gases or liquids. For 
example, it is well established for humid air accounting 
for which is important for realistic meteorological prog- 
noses, for salt water, hence oceanography, or even for flu- 
ids near the critical point met in engineering application. 
The paper also shows that the proposed generalization 
respects the conservation of the total mechanical energy. 
Obeying the conservation laws is of principle importance 
in a construction of consistent models, and a violation of 
this rule may lead to unphysical predictions. For strat- 
ified turbulent boundary layers this requirement lays in 
an explicit consideration of not only the mean profiles, 
but also of all relevant second-order, one-point, simulta- 
neous correlation functions of all fluctuating fields. 
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